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Finding the sigma pole by analytic extrapolation of 

7T7T scattering data 

Irinel Caprini, 

National Institute of Physics and Nuclear Engineering 
POB MG 6, Bucharest, R-077125 Romania 

be" 

O 
O 
(N 

u . 

Qh' We investigate the determination of the a pole from inr scattering data below 

the KK threshold, including the new precise results obtained from K e ^ decay by 
CS| ' NA48/2 Collaboration. We discuss also the experimental status of the threshold 

parameters a(j and 6{j and the phase shift <5q. In order to reduce the theoretical bias, 
we use a large class of analytic parametrizations of the isoscalar S-wave, based on 
expansions in powers of conformal variables. The a pole obtained with this method 
is consistent with the prediction based on ChPT and Roy equations. However, the 
theoretical uncertainties are now larger, reflecting the sensitivity of the pole position 
to the specific parametrizations valid in the physical region. We conclude that Roy 
equations offer the most precise method for the determination of the a pole from 
7T7T elastic scattering. 

> 

PACS: 13.75.Lb, 14.40.Cs 



1 Introduction 



The determination of the pole associated to the a resonance (or /o(600)) is known to be 
a difficult problem. The pole is situated deep in the complex plane, its influence in the 
physical region is masked to a certain extent by the nearby Adler zero and, until recently, 
the experimental data on mr scattering at low energies were quite poor. This explains why 
the values reported by PDG [1] for the mass and width of a cover a very large interval. 

During the last years, Chiral Perturbation Theory (ChPT) and Roy equations led to 
an accurate description of 7T7t scatering at low energies [2J [3] . In particular, the scattering 
length Oq and the effective range parameter b® of the isoscalar S"-wave given in [3]: 

a° = 0.220 ± 0.005, ^ = 0.276 d= 0.006, (1) 

have remarquably small uncertainties. 

The formalism based on Roy equations was shown recently [4] to control also the 
analytic extrapolation of the 7T7T amplitude in the complex plane, leading to precise values 
for the mass and width of o: 

M a = 441 +l 6 MeV, Y a /2 = 272 t\ 2 5 MeV. (2) 

In the standard method of detecting resonances, the experimental data on the partial 
wave with the quantum numbers of the resonance play an important role. Unlike this, 



the prediction (J2J) was obtained without using experimental data on the isoscalar S-wave 
at low energies: the amplitude was calculated below 800 MeV, and also in the complex 
plane, from Roy equations, using experimental input at higher energies and theoretical 
results on the pion-pion scattering [3J. Roy equations provide a very suitable framework 
in this case, compensating the lack of experimental data on im scattering at low energies 
by theoretical information. 

Recently [5], NA48/2 Collaboration measured the phase shift difference <5q — 5\ at low 
energies from K e ± decay, with a precision much greater than that of the older experiments 
P, [7j . This revived the interest in the determination of the scattering length a® and the 
pole associated to a by direct analytic extrapolation of the im scattering data. In j8] 
the authors propose a representation of the isoscalar S-wave £q(s) based on an expansion 
in powers of a conformal mapping variable. To account for the theoretical uncertainties 
related to analytic continuation, two parametrizations were considered, the difference 
between them being interpreted as a systematic theoretical uncertainty of the method. In 
the framework discussed in [8], the mass and width of a are obtained with an accuracy 
comparable to that quoted in (j2|). 

In the present work we focus on the problem of systematic uncertainties within this 
approach. We note that the class of functions used in [8], although based on a convergent 
expansion, is still quite narrow when the expansion is restricted to a few terms. By 
enlarging the class of admissible analytic functions used for fitting the data, the theoretical 
bias is reduced and a more realistic estimate of the uncertainties in the position of the a 
pole is obtained. In the present work we apply this idea, by using a large sample of analytic 
parametrizations of the tttt amplitude, suitable at low energies. A short description of the 
method and some results were given already in [9]. 

In the next section we discuss several parametrizations of the 7T7t isoscalar S'-wave, 
which satisfy analyticity and elastic unitarity. In the next two sections we apply these 
parametrizations for fitting the data on the phase shift 5®: in section 3 we consider only 
the data from K ei decay, and in section 4 we include data up to the KK threshold. From 
the admissible parametrizations of the isoscalar S'-wave we find the threshold parameters 
o° and &{] an d the location of the a pole. Our conclusions are summarized in section 5. 

2 The isoscalar S'-wave at low energy 

We consider the irir isoscalar S-wave £q(s), which is an analytic function in the s-plane 
cut along s > 4M% and s < 0. We assume that t[j(s) is the pure strong amplitude, where 
all the isospin breaking corrections are neglected. As in [2] -[3], we take for M^ and Mk 
the masses of the charged pion and charged kaon, respectively. 

Neglecting the inelasticity due to the 4ir channel below 1 GeV, unitarity implies that 
the relation 

Im -777 — 

is valid up to the threshold for KK production, s = 4M|-. From (jHJ) it follows that the 



ip{s), p(s) = Jl-4MZ/s, (3) 



function ip(s), defined by 



ip(s) -ip{s) 



is real in the elastic region: 

haip(s + it) = 0, AM 2 < s < am k, ( 5 ) 

and is related to the phase shift Sq(s) by: 

i>(s) = p(s)cot6° (s). (6) 

Since the amplitudes are analytic functions of real type, Eq. §5$) means that ip(s) has 
no discontinuity across the elastic unitarity cut. The definition (j4]) shows also that ip(s) 
has poles at the points where t®(s) has zeros. The amplitude is expected to vanish below 
threshold at a point sa, related to the so-called Adler zeros. ChPT to lowest order predicts 
sa = M 2 /2. Assuming that t®(s) does not have other zeros in the complex plane, the 
product (s — sa)iP{s) is analytic in the s-plane cut for s < and s > AM\. The effective 
range expansion amounts to expanding the function ip{s) in powers of q 2 = (s/4 — M 2 ) 
around the threshold q 2 = 0, where it is regular. However, the branch point s = limits 
the convergence of this expansion to the circle | q 2 \ < M 2 . 

2.1 Method of conformal mappings 

The domain of convergence of a power series can be enlarged by expanding the function 
in powers of a variable which conformally maps a part of the holomorphy domain onto the 
interior of a disk. The use of conformal mappings in particle physics was first discussed 
in [TOl ITT] ; in the context of the effective range expansion for partial waves a conformal 
mapping was used in |12j : more recently, the method was applied for the description of 
exclusive semileptonic B decays [13] and in perturbative QCD [14J. As shown in [TO] , 
the asymptotic rate of convergence of the series in the physical region is optimal if the 
amplitude is expanded in powers of the variable which mapps the entire holomorphy 
domain onto a disk. Since the disk is the natural convergence domain of the power series, 
the new expansion will converge in the whole analyticity domain, up to its boundary. 
Consider the variable 



W ( S) a ) = v (7) 

v/i + aJAM 2 K - s 



where a > is arbitrary. The function w(s,a) transforms the s-plane cut along s < 
and s > AM\ onto the unit disk |iu| < 1 in the complex plane w = w(s,a), such that 
u>(4M|-, a) = 1 and w(0, a) = — 1. In [8] the authors adopt the expansion 



M 2 



sa 



' B + Biw(s,a) + B 2 w(s 1 a) 2 + . 



M v ^~s 



(8) 



with the particular choice a — 1. In Eq. (jBJ), the first term in parantheses, added to 
the expansion in powers of w(s,a), compensates the singularity of p(s) at s = in the 



denominator of pi) , removing an unphysical singularity of £{j(s) on the real axis which 
would appear otherwise. 

A slightly different form was also used in [8] : 



i>(s) 






2s A 

M n 



B + B lW (s,l) + B 2 w(s,l) 2 + 



(9) 



where the factor (//jj — s) displays explicitly the energy where the phase shift 5® passes 
through 7r/2, according to (jSJ). This factor is useful for fitting narrow resonances but, as 
we shall show, it is not suitable for broad resonances like a. 

The power expansions in (jSj) and © converge in the disk \w\ < 1 and, for a large 
number of terms, these parametrizations are equivalent. However, when the series are 
truncated at a finite number of terms, (jSJ) and (jSj) lead to different results. This difference 
is interpreted in [8] as a systematic uncertainty of theoretical nature, which should be 
added to the statistical errors. In the present work we develop this idea, presenting other 
admissible analytic parametrizations of the amplitude. 

A first generalization is to expand ip(s) in powers of w — w(s, a), for an arbitrary a, 
as in ([8]). By varying a, one changes the point mapped to the origin of the u>-plane and 
the position of the intervals where experimental data are available. Some examples are 
shown in Fig. 1. As we shall see, the flexibility offered by the parameter a allows us to 
describe the peculiar structure of the isoscalar 5-wave near the inelastic KK threshold. 
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Figure 1: The disk \w\ < 1 in the complex plane w = w(s, a) defined in (|7j), for a = 0.36 
(left), a = 1 (center) and a = 4 (right). The thick segments indicate the regions where 
experimental data are available from K e4: decay [5]-[7j and the process irN -^ ttttN (cf. 
the compilation of data made in [IS]), respectively; the circle shows the a pole on the 
second Riemann sheet from [3]. 



2.2 Alternative procedure for ghost removal 

The singularity at s = of the phase space factor p(s) in (J3J) can be alternatively elimi- 
nated if the term ip(s) is replaced by a function which is analytic in the s-plane cut along 
s > 4M% and has the imaginary part equal to p(s) on the upper edge of the cut. In the 
context of effective range approximation, Chew and Mandelstam [16j defined such a func- 
tion, vanishing at threshold, by a once subtracted dispersion relation. For convenience, 



we consider the loop function of ChPT, J(s, M%), written as 



J( S , M^ = - + 



2 , P(s) 



In 



7T 7T 



l + p(s) 



which vanishes at the origin, J(0, M%) = 0, and satisfies the relation 

lm,J(s + ie,M 2 ) = p(s), s > AM 2 . 

If we define the function ipi(s) by: 

1 



t°o(s) 



ih{s)-J(s,M*y 



(10) 



;n) 



(12) 



the unitarity relation ([3]) and Eq. (fTTT) show that ipi(s) is real for AM 2 < s < 4M|-, where 
it is related to the phase shift <5q by 



Ms) = P(s) cot«Jg(s) + Re J(s, Af*) . 



(13) 



The reality property implies also that ipi(s) is analytic in the s-plane cut for s < and 
s > 4M|-, except for the pole at s = sa, and can be expanded as 



Ms) 



Ml 



s- S A 



B + B 1 w(s, a) + B 2 w(s, a) 2 + . . . 



(14) 



in powers of the variable (J7J). We remark that the compensating term 2saJM 2 ^[s ap- 
pearing in (jSJ) is no longer necessary in (}T4"j) . since the function J(s,M 2 ) is by definition 
regular at s = 0. 



2.3 .S-matrix factorization 

Other parametrizations of £q(s) are obtained by including some information about its 
behaviour near the KK threshold. We do this by expressing the S'-matrix element 



as a product 



S° (s) = 1 + 2ip(s)t° (s) 



Sq(s) = S Icst (s)S fo (s), 



(15) 
(16) 



where each factor satisfies elastic unitarity (|5' rest (s)| = |5/ (s)| = 1) below the KK 
threshold. The multiplication of the two ^-matrices amounts to the following addition 
rule for the corresponding amplitudes: 



where 



t 



to(s) = t rest (s) + tf (s) + 2ip(s) t rest (s) tf (s), 
Srcst(s) - 1 , , , Sf Q (s)-l 



(17) 



2ip(s) 



tfo(s) 



2ip(s) 



Crossing symmetry implemented by Roy equations [3] implies that the expansion of the 
partial wave amplitude around s = starts with 



t° {s) = t + t 1 s + t 2 s 3/2 + O(s' 



(19) 



where to is nonzero. In order to cancel the singularity of the factor p(s) at s = in 
(TT71) . either t rest (s) or t/ (s) must vanish at s — 0, but not both (since, cf. (IT91 . the 
full amplitude does not have a zero there). We choose to set t/ (0) = 0, taking for this 
amplitude the expression 



tfo( s ) 



kis 



K-s-hs J(s, Ml) - (k 2 + hs)J{s, M| 



2 ^' 



(20) 



where J(s, M 2 ) is defined in ( TTUi) and J(s, M^) is obtained replacing M w in p(s) by M^. 
We note that by taking 



1.01, fci = 0.08, fc 2 = -1.09, fc 3 = 1.16, 



(21) 



the modulus of the corresponding 5-matrix, Sf (s), is close, in the range 1M K < y/s < 
1.16 GeV, to the elasticity T]q(s) measured in [T7j, while for 



K 



-- 1.15, fci 
1.41, 



0.11, h 



0.39, 



0.03, 

= 1.72, 



:22i 



k x = 0.24, k 2 = -0.73, k 3 = 1.72, (23) 

|S/ (s)| is close, in the same range, to the upper/lower edges of the band of the elasticity 
?7g extracted from the decay J/ip -^ (fin-rc [18J. 

For our purpose, the specific form adopted for tf (s) is not a limitation, since the total 
amplitude contains the additional term t res t(s). Elastic unitarity, |S , res t('S)| = 1, implies 
that t rcst (s) can be written, for instance, as 



tr 



^rest(s) - ip{s) ' 



(24) 



with ip rest analytic in the s-plane cut along s < and s > 4M^, except for a pole at 
s — Si, where t res t( s i) = (from (|T7j) and (1201) it follows that si is close to the Adler zero 
sa)- Therefore, we can write V'rest(s) as: 



^rest(s) 



Ml 



S — S\ 



2s 



M n ^ 



— + B + B]w{s, a) + B 2 w(s, a) 2 + . . . 



(25) 



where w(s, a) is defined in ((7j). Alternatively, we can use for £ res t(<s) an expression similar 
to ( fT2l) . involving the function J(s,M 2 ). 

Other admissible parametrizations are obtained if we assume that t rest (s) is almost 
regular near s = 4M|-. Since the next branch point, at s — 4M^, is known to have a weak 
effect, we can neglect at low energies the right hand cut of ip TCS t(s), and expand it as 



V'r 



Ml 



Sl 



2si 

M^yfs 



+ B + B\W\{s, a) + B 2 wi(s, a) 2 + . . . 



(26) 



where the variable 

w 1 {s,a)= V f~ a , a>0, (27) 

maps the s-plane cut only for s < onto the disk \wi\ < 1 of the complex plane w\ = 
Wi(s, a). 

The expressions given in this subsection are examples of possible analytic parametriza- 
tions of the amplitude at low energies. Following an idea of Dalitz and Tuan [T9], in 
some phenomenological analyses [20] the individual 5*-matrices in the product f|T6|) are 
associated to specific resonances. However, in our work we use the factorization ([TBI 
only for mathematical purposes: by isolating a factor with a rapid variation near the 
KK threshold, we expect a better convergence for the expansion of the remaining part, 
treat (s)- This part, which is fixed by the low energy data, contributes to both the elasticity 
r)o(s) = |5o(s)| and the phase shift 5$ above the KK threshold. So, the behaviour of tg 
above this point is left free in our fits. 

3 Fits of the data from K e ^ decay 

We consider first the data on the difference 5q—5\ measured below 0.4 GeV from K^ decay 
[5]- [TJ. The P-wave phase shift 5\ is known with precision in this energy range [3l [15], 
allowing an accurate extraction of 5$. As shown recently [21], the phase shift measured in 
K e4: decay differs from the pure strong phase shift 5q(s) by an isospin correction overlooked 
so far, accounting for the differences between the masses of the charged and neutral 
mesons, and between the quark masses m u and mj. The correction evaluated in ChPT 
to one- loop reads [21] : 



where p(s) is defined in (J3J) and 



^ = Ml-Ml 0l po(s) = Jl-AM^/s, R=— , m = (m u + m d )/2. (29) 

With the estimate R=37 ± 4 given in [21], the correction A[5q(s)] amounts to a fraction 
of a degree in the whole experimental range. This correction was subtracted from the 
phase shift derived from K e4 data, in order to obtain the pure strong phase shift Sq(s). 

The total number of points from the K ei experiments is 21 (5 points from [6], 6 from 
[7J and 10 from [5]). As in [8], we increased the experimental error on the last point in 
[7J by 50%. For the 10 data from the NA48/2 experiment we used the covariance matrix 
published recently in [5]. We fitted these data with the parametrizations described in 
section 2. 

In our analysis, the positive number a, specifying the conformal variables and 
f[27j) . together with the parameters k and ki appearing in ([201 . represent the input which 
defines an admissible class. To account for the uncertainty in the position of the Adler 
zero, sa was varied between 0.4M^ and 0.6M^. In each admissible class, the coefficients 
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Figure 2: Left: phase shift 8® derived from K e ^ decay, fitted with the 16 parametrizations 
described in the text. Right: extrapolation of the fits above the experimental range. 



Bi of the expansion in powers of the conformal variable are free. They are determined by 
fitting the low energy data. 

We investigated a large class of combination of input parameters, from which we 
retained 16 admissible parametrizations: the first three are based on Eqs. (J3J) and (jSJ), 
with a = 1, a = 0.36 and a = 4, respectively. The value a = 0.36 is special, since, as seen 
in Fig. 1, it maps the experimental range relevant in K ei decay onto a symmetric interval 
around the origin w = 0. According to general theorems |12j . this variable gives the best 
approximation of data in the experimental range. On the other hand, for a greater than 
1 (for instance, a = 4) the variable w(s,a) maps the region close to the branch point 
s = 4M|- near the origin of the w-plane (see right panel of Fig. 1). As we shall see, the 
expansion ([8]) is then more suitable at higher energies. 

The next three parametrizations are based on Eqs. (TT2l) and (TT4l) . with the same 
choices a = 1, a = 0.36 and a = 4. The difference from the previous fits is the way of 
eliminating the singularity at s = of the phase space p(s) in the denominator of t^(s). 

The remaining ten parametrizations are based on the S'-matrix factorization discusssed 
in subsection 12.31 In three of them we take the parameters k and ki from fT22|) and the 
expansion fT26|) . with a = 0.36, a = 1 and a = 4, and in the next two we use the same 
parameters k and ki and the expansion fT2oT) . with a = 0.36 and a = 1. In other two 
cases we use the parameters n and fcj from (T2"3"|) and the expansion (12"6"|) with a = 0.36 
and a = 1, while in the last three parametrizations we take the parameters k and fc, from 
( 12T|) . using either the expansion (|25|) with a = 1 and a = 0.36, or the expansion (1261) with 
a = 0.36. 

We obtained good fits of the 21 experimental points with 2 free parameters, Bq and 
Bi, in the expansion in powers of the conformal variables. The values of x 2 are very 
similar for all the fits, although the parametrizations are quite different. The values of 
X 2 and the optimal parameters are given in Table 1, for sa (or s\) fixed at 0.5M^ (note 
that if Si = 0.5M^, the position sa of the Adler zero, resulting from the fits, is slightly 
different: for the last 10 fits given in Table 1, sa varied between 0.42M^ and 0.47M^). 



Nr. 



r 



Bn 



Bx 



10 

u o 



/g7 (MeV) 



1. 


21.7 


7.5 


-15.1 


0.216 


0.278 


459 + 259 i 


2. 


21.5 


14.6 


-12.4 


0.214 


0.282 


445 + 259 i 


3. 


21.9 


-16.4 


-37.1 


0.217 


0.275 


473 + 261 i 


4. 


20.9 


7.7 


-20.2 


0.212 


0.287 


412 + 237 i 


5. 


20.6 


17.2 


-16.5 


0.210 


0.292 


401 + 231 i 


6. 


21.2 


-35.0 


-60.3 


0.214 


0.284 


422 + 246 i 


7. 


21.5 


14.6 


-14.8 


0.214 


0.281 


443 + 262 i 


8. 


21.7 


5.8 


-18.8 


0.215 


0.278 


455 + 261 i 


9. 


21.8 


-25.3 


-47.9 


0.216 


0.276 


465 + 260 i 


10. 


21.6 


15.2 


-12.8 


0.215 


0.280 


451 + 264 i 


11. 


21.8 


7.8 


-15.5 


0.216 


0.277 


466 + 264 i 


12. 


21.5 


15.0 


-15.3 


0.214 


0.281 


448+ 264 i 


13. 


21.7 


5.9 


-19.5 


0.216 


0.277 


459 + 262 i 


14. 


21.8 


7.8 


-15.5 


0.216 


0.277 


465 + 263 i 


15. 


21.6 


15.1 


-12.8 


0.215 


0.280 


450 + 263 i 


16. 


21.4 


14.5 


-14.8 


0.214 


0.282 


443 + 261 i 



Table 1: Results of the fits of the data from K e4 decay [5]- [7J, using the 16 parametrizations 
described in the text. 



The values of x 2 decrease by about 0.4 if we take into account the theoretical uncertainty 
associated to the isospin correction [21 J. For simplicity we indicate only the central values 
of the parameters, omitting the statistical errors. 

The quality of the fits is seen in Fig. 2. Although the fits are almost indistinguishable 
in the experimental range, they exhibit large differences when extrapolated to higher en- 
ergies. This illustrates the well-known phenomenon of instability of analytic extrapolation 
[22] . We note that the lowest curves in the right panel of Fig. 2 correspond to the fits 4, 
5 and 6 in Table 1, obtained with the parametrization ffT2l - ffT4l) , In particular, the fits 
no. 4 and 5, corresponding to the choices a = 1 and a = 0.36 in ()14[) . exhibit a plateau 
at low values of $j(s). The increase of the phase shift required by the high energy data 
(see below) is obtained, for instance, with the choice a = 4 in the expansion (jSJ), or by 
using parametrizations based on the S'-matrix factorization described in subsection 2.3. 

In Table 1 we give for each fit the central values of a{] and 6° and the position s a of the 
a pole on the second Riemann sheet, obtained by analytic extrapolation to the threshold 
and into the complex plane (as shown in [I], s a is the solution of the equation S(s a ) = 
on the first sheet). 

Taking the average of the 16 admissible values for m a = J~s^. = M ff —iT ff /2, we obtain: 



M„ 



447 ± 7 (stat) j^jj (syst) MeV, T a /2 



258 ± 6 (stat) t*° (syst) MeV, 



(30) 



where the systematic error is defined, as in [8] , such as to cover all the admissible fits (the 
uncertainty in the Adler zero produces a small error, of about 4 MeV in M a and 3 MeV 
in IV/2). 
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Figure 3: Positions of the a pole obtained by the analytic extrapolation of the parametriza- 
tions used for fitting various sets of data, compared with Refs. [1], [S] and [23] (from the 
last reference we show the value obtained with isospin corrected K e ± data). 



The pole positions given in Table 1 are shown in Fig. [3j together with the results 
reported in [1] , [8] and [23] . Note that the three isolated points, with small values of M a 
and r ff , correspond to the fits 4, 5 and 6 in Table 1, which are based on the parametrization 
(TT2l-(fT4"l). As we mentioned, they lead to a bad behaviour of the phase shift at higher 
energies, in spite of the fact that they provide very good fits of the data from K e ± decay. 
From Table 1 it is seen that these parametrizations (especially 4 and 5) give low values 
for Oq and large values for &{]. 

In this Section we obtained good fits of the phase shift measured from K ei decay. 
However, the extrapolation of the phase shift above the experimental region is not ac- 
ceptable for many of them. Therefore, we can not take the results of this Section as final. 
In particular, the narrow range spanned by most of the widths IV in Fig. [3] may signal a 
bias. In the next Section we shall improve the description of £{j(s) by including data on 
the phase shift at higher energies. 



4 Inclusion of high energy data 

The difference <5q(s) — Sq(s) is measured at s — M\ from the decay K — > im [21]. However, 



in this case the radiative corrections are very large [25] and the extraction of the strong 
phase shift 5$ is still uncertain [26]. For this reason, we shall not use as input in our 
analysis this datum. 

Experimental data at higher energies are available from the irN — > ttttN process 
[27] [25] . We considered two sets of data below the KK threshold: 

• set I, which consists from 40 data points: 21 from K ei decay [S]-[S] and 19 from the 



10 



Nr. 



r 



XKe4 



B n 



B l 



By 



b° 



/^(MeV) 



1. 


37.7 


24.3 


7.8 


-23.5 


-20.6 


0.233 


0.261 


486 + 312 i 


2. 


32.9 


22.8 


-20.3 


-61.5 


-23.9 


0.226 


0.271 


462 + 298 i 


3. 


32.6 


22.6 


-37.3 


-84.6 


-29.6 


0.225 


0.272 


461 + 296 i 


4. 


32.7 


22.2 


1.3 


-52.3 


-39.9 


0.222 


0.271 


458 + 265 i 


5. 


33.9 


21.8 


-52.8 


-105.8 


-26.0 


0.207 


0.283 


442 + 234 i 


6. 


32.4 


22.1 


9.3 


-12.0 


- 


0.220 


0.274 


457 + 281 i 


7. 


33.7 


22.9 


-2.2 


-20.3 


- 


0.228 


0.272 


454 + 303 i 


8. 


33.9 


21.8 


13.1 


-12.2 


- 


0.215 


0.278 


454 + 267 i 


9. 


38.0 


21.8 


5.2 


-19.9 


- 


0.213 


0.278 


456 + 256 i 


10. 


32.1 


21.5 


18.8 


-12.9 


- 


0.216 


0.281 


441 + 272 i 


11. 


31.9 


21.6 


12.4 


-13.9 


- 


0.216 


0.279 


446 + 267 i 


12. 


32.1 


21.9 


7.4 


-15.5 


- 


0.219 


0.276 


450 + 274 i 


13. 


32.7 


22.6 


-9.5 


-29.3 


- 


0.225 


0.272 


455 + 295 i 



Table 2: Results of the fits of the 40 data points of set I, using the 13 parametrizations 
described in the text. 



CERN-Munich experiment |27j ; 

set II, which consists from 32 data points: 21 from K ei decay [S]-[0] and a collection 
of 11 data points from ttN —> ttttN [2B], given in Eq. (2.13) of 



As in the previous section, we investigated a large number of parametrizations, but 
rejected many of them since they gave bad fits. For instance, the choice a = 0.36 (or other 
values a < 1) in the expansions (jHD and (TT4"j) was not admissible, leading to high values of 
X 2 (such parametrizations can not exhibit the rapid increase of the phase shift above 900 
MeV). Also, the parametrization (Q, considered in [8], proved to be not acceptable: with 
3 free parameters, /io, B and Bi, we obtained \ 2 — 45.6 for the 40 points of set I, and 
X 2 = 34.6 for the 32 points of set II. We recall that expressions which display the energy 
where the phase shift passes through n/2 are often used for fitting narrow resonances. 
However, they are not suitable for broad resonances like a. 

We finally retained 13 admissible parametrizations: the first three are based on Eqs. 
dlj and ([8]), with a = 1, a = 4 and a = 6, respectively. The next two are based on Eqs. 
( fl2|) and (fl4|) . with the choices a = 1 and a = 4. The remaining eight parametrizations 
are based on the S'-matrix factorization discusssed in subsection 12.31 in four cases we use 
the values k and k t from (12 ip . and either the expansion (1251) with a = 1, a = 4 and 
a = 0.5, respectively, or the expansion (126]) with a = 1. In the next three cases we use 
the parameters k and k{ from (1221) . and either the expansion (125]) with a = 0.2, or the 
expansion (126]) . with a = 0.5 and a = 1. Finally, in the last parametrization we take the 
values of k and k,i from (T2"3"]) and the expansion (T261) with a = 4. 

For the first five parametrizations we used 3 nonzero coefficients, B , B\ and B 2 , in 
the expansion in powers of the conformal variables, and for the last eight we obtained 
good fits with 2 nonzero coefficients, Bo and B\. As in the previous Section, we took into 
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Nr. 



r 



XKe4 



B n 



B l 



By 



b° 



^(MeV) 



1. 


28.6 


22.8 


5.4 


-29.0 


-20.9 


0.225 


0.267 


494 + 279 i 


2. 


24.8 


22.4 


-26.6 


-71.8 


-27.1 


0.223 


0.272 


470 + 287 i 


3. 


24.4 


22.3 


-45.9 


-97.4 


-33.0 


0.223 


0.273 


469 + 285 i 


4. 


25.1 


22.0 


-0.3 


-56.6 


-41.3 


0.218 


0.273 


466 + 251 i 


5. 


24.6 


22.2 


-56.6 


-111.1 


-26.9 


0.205 


0.284 


446 + 229 i 


6. 


23.9 


21.8 


7.4 


-16.3 


- 


0.215 


0.277 


466 + 259 i 


7. 


24.3 


22.5 


-6.2 


-25.3 


- 


0.224 


0.272 


462 + 293 i 


8. 


24.7 


22.1 


12.5 


-16.9 


- 


0.208 


0.281 


463 + 242 i 


9. 


26.6 


22.9 


1.0 


-28.8 


- 


0.204 


0.282 


465 + 225 i 


10. 


23.2 


21.6 


20.1 


-16.8 


- 


0.210 


0.284 


448 + 254 i 


11. 


23.2 


21.8 


11.6 


-18.7 


- 


0.209 


0.283 


453 + 246 i 


12. 


23.1 


21.8 


4.8 


-20.8 


- 


0.213 


0.279 


458 + 253 i 


13. 


23.5 


22.0 


-21.1 


-43.3 


- 


0.219 


0.274 


466 + 272 i 



Table 3: Results of the fits of the 32 data points of set II, using the 13 parametrizations 
described in the text. 



account the uncertainty in the position of the Adler zero, by allowing in each case sa to 
vary between 0.4 M 2 . and 0.6 M%. 

The results of the fits (for sa or s\ fixed at 0.5M 2 ) are presented in Tables 2 and 3, for 
the sets I and II, respectively. For completeness, we show also in column 3 the contribution 
to the x 2 °f the 21 points from K e ^ decay. By comparing Tables 2 and 3 with Table 1 we 
notice that the description of the K e ± data, measured by their contribution to the total x 2 , 
is now slightly worse than in the fits restricted to the data from K e 4 decay: constraining 
the behaviour at high energies leads to a small deterioration of the description of the low 
energy data. The quality of the fits is shown in Fig. HJ and in Fig. 5 we show an expanded 
view of the energy region covered by K e ± decay. As expected from the previous discussion, 
the various parametrizations are no longer indistinguishable, as were those obtained from 
fitting only the K e 4 data in Fig. 2. 

Taking the average of the values of the threshold parameters a[] and 6q given in Tables 
2 and 3, weighted with the ratio xV^dof, we obtain, for the two sets: 



0.220 ± 0.005 (stat) ± 0.013 (syst) ± 0.003(s A ) (I) 
0.215 ± 0.005 (stat) ± 0.011 (syst) ± 0.003(s A ) (II) 



(31) 



and 



b° = 0.275 ± 0.006(stat)+J;[| ( ] ) 4(syst) ± 0.004 (s A ) (I) 
b° = 0.277 ± 0.006(stat)+J;J°J(syst) ± 0.004 (s A ) (II). 

If we combine these determinations we obtain the values 



(32) 



0.218 ±0.014, 



b° 



0.276 ±0.013, 



(33) 
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Figure 4: Left: fits of the data in set I (K eA data [5]-[7j plus the CERN- Munich data [27] 
below the KK threshold). Right: fits of the data in set II {K eA data plus a selection of 
data from ttN — » mtN, given in |15j). 



shift 5%(s 



which are fully consistent with the results obtained from ChPT and low energy theorems 
for 7ttt scattering [3], quoted in ([T]). 

Before discussing the results for the a pole, let us make a few comments on the phase 
shown in Fig. 4. It was advocated in some papers, for instance [23], that the 
phase shift of the isoscalar S'-wave exhibits a "hump" at energies around 800 MeV, before 
starting the rapid increase near the KK threshold. As seen from Fig. 4, no hump is seen 
in the fits of the data in the set I, while a weak hump appears only in a few fits of the data 
in set II. This proves that the hump seems to be an artefact of special parametrizations 
used for fitting the data (in particular, the expression (jUj) displays a pronounced hump 
if a small number of coefficients Bi is kept in the expansion). As discussed in [29J, this 
shape is in conflict with the forward dispersion relation for the tctt amplitude of isospin 
1 = 0. 

It is of interest to calculate with our parametrizations the value of Sq at \fs = Mr. We 
recall that the phase shift difference 5q(s) — 5q(s) at this energy can be extracted from the 
decay K — > titt. However, as we mentioned in section 4, in this case the isospin breaking 
corrections are large and the extraction of the strong phase shift 5$ is still unclear. For 
this reason we did not use this information as input in our fits. Using the parameters 
given in Tables 2 and 3 and taking the averages of the admissible values in the two sets 
we obtain: 






1.7° 



38.9° ± 0.6° (stat)I^o (syst) (I) 



2.6° 



40.4° ± 0.8° (stat)li;g (syst) (II). 



(34) 



ID 



48.7 ± 4.9°, which is 



In [S] the authors used as input in their fits the value <5q(M| 
significantly larger than the output values given in (134"]) . 

The phase shift at 0.8 GeV is also of interest, since it is a key input in solving Roy 
equations |2J. We recall that in [2] the range is r5°(0.8 GeV) = 82.3 ± 3.4° was adopted as 
input, while in [3] the more conservative choice <5g(0.8GeV) = 82.3l 4 o was made. From 
the parametrizations discussed above and the parameters Bi given in Tables 2 and 3, we 



13 



30 r 



25- 



20- 



15 



10- 



Geneva-Saclay (1977) 

E865(2001) 

NA48 (2008) 



5- 






# 



m 



*& 



.1- 



>-" 



"f 



..*? 



.,.<- 



30 

25 

20 
) 

) 

15 

10 
5 




Geneva-Saclay (1977) 

E865(2001) 

NA48 (2008) 



** 



+ 



,.,,i 



aj 



<C 



J%* 



,.•*• 



,i- 



0.3 



0.33 0.36 

Energy (GeV) 



0.39 



0.3 



0.33 0.36 

Energy (GeV) 



0.39 



Figure 5: Expanded view of the low energy region from Fig. HI 



obtained the average values: 



$(0.8 GeV) = 81.8° ± 0.6° (stat) ± 1.3° (syst) 



+3.3° 



$(0.8 GeV) = 85.9° ± 0.7° (stat)^ (syst) 



(I) 
(II), 



(35) 



which are consistent with the ranges adopted in [2] and [3]. 

We turn now to the predictions for the a pole obtained by the analytic continuation of 
the parametrizations considered above. The a pole positions given in Tables 2 and 3 are 
shown in Fig. [3J The three isolated points obtained with some fits of the K ei data are no 
longer allowed, but in the same time the tight correlation exhibited by the other fits of 
the K e 4 data is now softened. This is due to the fact that the description of the K ei data, 
measured by their contribution to the total x 2 5 is slightly worse than in the fits restricted 
only to the K e 4 data, and the various parametrizations are not as indistinguishable at low 
energies as in Fig. 2. 

The averages of the values given in Tables 2 and 3, respectively, weighted with the 
corresponding x 2 /N Aoi) give: 



M„ 



278 ± 6(stat)±2(syst) MeV (I) 



455±6(stat)+g(syst)MeV, T a /2 

463±6(stat)+^(syst)MeV, r CT /2 = 259 ± 6 (stat) t^ (syst) MeV (II). (36) 



As above, the systematic errors cover the values in the admissible samples. The uncer- 
tainty in the position of the Adler zero sa has now a smaller effect, of 2 MeV for M a and 
1 MeV for V a j1. Alternatively, we can define the central values from the optimal fits with 
the lowest x 2 (fit no. 11 in Table 2 and fit no. 12 in Table 3). This procedure gives 



M„ 



1-40 



446 ± 6(stat)7X (syst) MeV, T a /2 
458±6(stat)+^(syst)MeV, T a /2 



267 ± 6(stat)+JJ(syst) MeV (I) 

253 ± 6(stat)+2g(syst) MeV (II). (37) 



However, this definition is not very sharp since, as seen from Tables 2 and 3, there are 
several fits with very close values of x 2 , which lead to different values for M a and T a . 



14 



O.N 



O.d 







0.4 



0.2- 



- ll. 1 - 




1 ' 1 ' 1 










, 


• Cohen (1980) 
■ Etkin (1980) 
BES (2005) 
Fits 



1.05 1.1 

Energy (GeV) 



1.15 



Figure 6: Elasticity 77g obtained with the 13 parametrizations of the data in set I (cf. 
Table 2), compared with the experimental data from [T7]-|18j. 



Eqs. (13T)]) and (I3"7j) represent our final results for the a pole position, obtained using 
the data on K e 4 decay and two, rather complementary, sets of scattering data at higher 
energies. The differences between them indicate the sensitivity of the pole location to the 
behaviour of the phase shift near the KK threshold. 

The comparison of ( |36|) and ( 1371) with (T5]) shows that the analytic extrapolation of 
experimental data leads to values for the mass and width of a which are consistent with 
ChPT and Roy equations, but have larger theoretical uncertainties. The errors are pro- 
duced by the well-known instability of analytic continuation [22]: functions very close 
along a limited part of the boundary may differ drastically outside the initial range. We 
illustrate this feature in Fig. 6, where we show the elasticity 77° calculated with the 13 
parametrizations fitting the data of set I, given in Table 2, extrapolated above the in- 
elastic threshold. Also, in Fig. 7, we show the real and the imaginary parts of the same 
parametrizations along a range covering a part of the left hand cut. Note that we are 
calculating now the amplitude on the cuts of the s-plane, while, strictly speaking, the 
expansions in powers of conformal variables converge only at points inside the analytic- 
ity domain [TU]. However, the expansions are truncated at low orders and are far from 
the asymptotic regime, so we may view them as effective parametrizations which have a 
meaning also on the boundary. 

We recall that information about t%{s) along the left hand cut (rigorously speaking, 
only on a part of it) is obtained from crossing symmetry. The dominant contribution is 
given by the p resonance, which does not make a narrow peak. Of course, in the present 
approach crossing symmetry is not explicitly implemented. Fig. 7 shows the differences 
in the values on the left hand cut of amplitudes which are almost indiscernable in the 
physical region. Since the a pole is rather close to the left hand cut, the spread in the 
positions of the a pole shown in Fig. [3] is not surprising. 
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Figure 7: Rei[j(s) (left) and Ihi£q(s) (right) obtained by the extrapolation of the 13 
parametrizations of the data in set I. 



5 Conclusions 

The new accurate data on ttit scattering at low energies obtained from K e ^ decay by the 
NA48/2 Collaboration [5] revived the interest in finding the a pole by the standard method 
used for narrow resonances, i. e. by the analytic extrapolation of a suitable parametrization 
of the partial wave with the quantum numbers of the resonance. 

In the present work we extended the investigation done in [8], by using a larger class 
of analytic functions for the parametrization of the titt isoscalar S'-wave at low energies. 
The purpose was to reduce the theoretical bias and to provide a more realistic estimate 
of the systematic uncertainties on the pole position. 

Our analysis shows that, in spite of the remarkable accuracy of the new data obtained 
from K e 4 decay [5], the inclusion of data at higher energies is necessary in order to 
reduce the theoretical bias and to exclude parametrizations which do not have a suitable 
behaviour above the experimental range. 

The values (155]) represent our prediction for the mass and width of a, obtained by the 
analytic extrapolation of a large number of admissible parametrizations of the isoscalar 
S'-wave. We present separately the results obtained by fitting with the same parametriza- 
tions the two sets of data from the process 7rN — > ttttN, in order to illustrate the sensitivity 
of the pole position to the behaviour near the KK threshold. We emphasize that, although 
we used a large number of parametrizations, the admissible sample is still limited. There- 
fore, the procedure is not entirely model independent (for a parametrization-free method 
for the detection of resonances from error-affected data given along a finite range see [30J). 

Our results (1361) are consistent with the mass and width of a obtained from ChPT and 
Roy equations, quoted in ([2]). However, the method employed here has larger theoretical 
uncertainties due to the phenomenon of instability of analytic extrapolation from a part 
on the boundary [22]: the differences between the various parametrizations are amplified 
by the extrapolation from the physical region to a distant point in the complex plane. 

The extrapolation error can be kept under control by using additional information 
about the physical amplitude, besides the low energy experimental data. In the method 
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based on Roy equations, this information is provided mainly by crossing symmetry and 
low energy theorems for tttt scattering [3]. As shown in [I], this tames the instability 
of the extrapolation to the a pole, leading to the small errors quoted in ([2]) (for a de- 
tailed discussion see also [21] )• We conclude that Roy equations provide the most precise 
determination of a from nn elastic scattering. 
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